Piezoresistive anisotropy of percolative granular metals 
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The piezoresistive response of granular metals under uniaxial strain is strongly dependent on the 
concentration of the conducting phase. Here we show that the piezoresistive anisotropy is reduced 
as the system approaches its percolation thresold, following a power law behavior in the critical 
region. We propose a simple relation between the conductance and the piezoresistive anisotropy 
which could be used in relation to real materials and notably to the thick film resistors. PACS 
numbers: 72.20.Fr, 72.60.+g, 72.80. Ng 
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Disordered systems constituted by random dispersions 
of metallic and insulating phases, also commonly referred 
to as granular metals, show transport properties strongly 
dependent of the concentration x of the metallic phase 
in the composite. For x larger than a critical concen- 
tration x c , transport is metallic-like while for x — > x c 
the system undergoes a metal-to-insulator transition. In 
the critical region (x — x c ) <C 1, the conductance G fol- 
lows a power-law behavior G ~ (x — x c ) where t is the 
critical exponent. The minimal model of a theoretical 
descriptions of transport in granular metals is based on 
extensions of ^percolation theory to random-resistors net- 
work models.™ Numerical simulations of such percolative 
resistors show that as x approaches x c from above, fluctu- 
ations of the microscopic (bond) currents get enhanced. 
It has been shown recently that this phenomenon can be 
spoiled to get more accurate experimental determinations 
of x c compared to common conductance measurements. □ 

The enhancement of bond current fluctuations is ba- 
sically also at the origin of the peculiar behavior of the 
piezoresistance effect RS X ► Xq. u Under an uniaxial im- 
posed strain, in fact, the piezoresistive response is highly 
anisotropic when x ~ 1, being maximum for external 
electric fields applied along the direction of the imposed 
strain (longitudinal response) and minimum along the or- 
thogonal directions (transverse response). This is due to 
the fact that for x ~ 1 the microscopic bond currents are 
basically all in the direction of the external electric field. 
Instead, lower values of x force the bond currents to have 
components also along directions perpendicular to the ex- 
ternal fields so that the anisotropy of the piezoresistive 
response is gradually reduced until at x = x c the longi- 
tudinal and transverse responses become equal.El Hence, 
at the percolation thresold the piezoresistive effect is per- 
fectly isotropic. 

The transport properties of granular metals under uni- 
axial strain could in principle help to characterize the 
vicinity to the percolation thresold with an higher accu- 
racy compared to that typical of common conductance 
measurements. The interest on the piezoresistive re- 
sponse of granular metals however concerns also more 
applicative problems. In fact, a particular class of gran- 



ular metals made of Ru02, IrC>2 or other metal-oxide 
granules embedded in a glassy matrix, also commonly 
known as thick-film resistors (TFRs), are characterized 
by a quite large piezoresistive response, which is usedj 
fabricating highly sensitive pressure and force sensors. 
It is thought that intergrain tunneling processes are re- 
sponsible for the TFRs high strain sensibility.ErQ A char- 
acteristic common to most of the TFRaJs. that their 
piezoresistive response is quasi-isotropic,&El'El indicating 
that the values and directions of the microscopic cur- 
rents are highly fluctuating functions and that the TFRs 
are quite close to the percolation thresold. This is con- 
firmed also by the observation o£, power-law dependence 
of G upon R11O2 concentrations.u 

In this communication, we compute the piezoresistive 
response of different random resistor network models as 
a function of bond or site concentrations. We show 
that close to the percolation thresold, the piezoresistive 
anisotropy follows an universal behavior governed by an 
exponential law. Our prediction could be easily experi- 
mentally tested, and could bring important informations 
on the topology of the percolative network and the criti- 
cal point of TFRs. 

To model the piezoresistive behavior of TFRs, or more 
generally of percolative granular metals, we consider a 
three-dimensional cubic random-resistor network whose 
bond conductances either zero or proportional to a tun- 
neling exponential factor exp(— 2dj/£), where di is the 
tunneling distance between two neighbouring sites along 
the i direction and £ the localization length. In our 
model, the network bond directions are parallel to the 
orthogonal x, y, z axes, and in the absence of applied 
strains the tunneling distances arc all equal: di = d, 
i = x,y, z. In this situation, the network is isotropic 
and the total conductance Gi = G depends only on the 
specific bond conductance distribution. To study the 
effect of an applied strain, let us assume that the re- 
sistor network is embedded in an homogeneous elastic 
medium and that the elastic coefficients of the network 
and the medium are equal. Hence, when a strain field 
£a (i — x, y, z) is applied to the sample, the conductance 
changes to Gi — G + SGi, where the variation SGi can be 
expressed in terms of the (intrinsic) conductivity varia- 
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FIG. 1. Monte Carlo calculations of longitudinal (Ty) and 
transverse (T±) piezoresistive coefficients of bond percola- 
tion (a) and site percolation (b) cubic network models. The 
calculations have been performed for networks of L 3 sites 
with L — 30 and L = 40 by numerically solving itera- 
tively the Kirchhoff node equations (for more details see Ref. 
fy. Filled symbols refer to bond conductances with tun- 
neling exponent 2d/£ = 4, while open symbols are the re- 
sults of uniform random distribution 2 < 2d/( t < 6. In (a), 
T|| and T± become equal at the critical bond concentration 
p = p c ~ 0.248. In (b), piezoresistive isotropy is reached at 
the critical site concentration x — x c ~ 0.314 (which corre- 
sponds to p c — x c ~ 0.098). 



tion 5ai and a geometric (extrinsic) factor:! 
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where a is the unstrained conductivity and the indexes 
i, j, k assume the values x, y, and z with cyclic permu- 
tations. In the above equation, the variation biii is due 
to the tunneling distance change d — > dj = d(l + en) 
which leads to exp(2<i/£) — > exp(2d/£)(l + 2deu/^) for a 
conducting bond along the i direction. 

To express Sai/a as a function of a general en, let us 
first imagine that an uniaxial strain along the x direction 
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0) is applied to the sample. In this 



situation, the conductivity along the x axis will in general 
be different from those along the y and the z axes which, 
by symmetry, are instead equal. Therefore, up to linear 
terms in e, a x = a — aT»e and a y = a z = a — aT±e where 
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are the longitudinal and transverse piezoresistive coeffi- 
cients, respectively. The above reasoning can be repeated 
for uniaxial strains along the y and the z axis and, since 
the problem is linear, for a general strain field ea the 
conductivity variations 5ai reduce to: 
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The above equation permits to express all the piezoresis- 
tive coefficients = —5ai/ejja or, by using Eq.(^j), 
the corresponding piezoresistive gauge factors K^j — 
—SGi/ejjG (that are the commonly measured quantities) 
in terms of r m and T±. For example, in a typical experi- 
ment with a thick cantilever beam having the main axis 
directed along the x direction, the strains are approxima- 
tively e xx — e, e yy — —ve, and e zz = —u's, where v and 
v' are the Poisson ratios of the cantilever and the resistive 
sample, respectively.E3 By using Eqs.(|l],|]), the longitudi- 
nal (-Kl) and transverse (Kt) piezoresistive gauge factors 
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K yx = -(i + r l{ )v- (i-r x )(i -v'). 
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Hence from a measurement of K xx and K yx it is possible 
to extract the intrinsic piezoresistive coefficients T» and 
T±. As we are going to show below, the knowledge of T|| 
and T± permits to characterize the nature of transport 
of percolative piezoresistors. 

As already pointed out in the introduction, the piezore- 
sistive response depends in a non trivial way on the vicin- 
ity of the resistor network to its percolation critical point. 
Consider in fact Fig. la in which we report numerical 
Monte Carlo calculations of Th and T± for bond perco- 
lating random resistor networks. In this model, a frac- 
tion p of bonds distributed at random has conductance 
exp(— 2d/£), while the remaining fraction 1—p has zero 
conductance. In the figure, filled symbols refer to a model 
in which the tunneling exponent 2d/£ is set equal to 4, 
while the open circles are obtained by treating £ aa a ran- 
dom variable in such a way that 2 < 2d/£ < 6.B When 
there are no missing bonds (p = 1) the piezoresistive 
response is anisotropic with the longitudinal coefficient 
much larger than the transverse one, reflecting the fact 
that the current flows mainly along paths parallel to the 
direction of the voltage drop. This is more evident for 
the 2d/£_ = 4 data (filled circles) for which T\\ = 2d/£ = 4 
and r_i_ = at p = 1. If fluctuations in the values of 2d/£ 
are allowed (empty circles) even at p = 1 the microscopic 
currents develop a component ortogonal to the voltage 
drop direction leading to rii < 2d/£ and T± > 0. When 
bonds are removed (p < 1), the longitudinal response gets 
reduced and at the same time the transverse piezoresis- 
tive coefficient is enhanced. This trend persists all the 
way down to the percolation thresold (p = p c ~ 0.2488) 
where Tu equals T± and the system shows perfect piezore- 
sistive isotropy. As shown in Fig. lb, the same qualita- 
tive behavior holds true also for site percolation models 
in which a concentration x of sites is removed at random 
with remaining conducting bonds having 2d/£ = 4 (filled 
squares) or 2 < 2d/£ < 6 (open squares). 

From the results of both bond percolation and site per- 
colation models of Fig.l we can infer that as a general 
rule 3T 1 ii — > r_i_ when the resistor network approaches its 
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FIG. 2. Piezoresistive anisotropy factor as a function of 
p — p c for the bond (circles) and site (squares) percolation 
data of Fig. 1. In the inset the same data of the main figure 
are plotted in a log-log scale. 



percolation thresold. This behavior is better analyzed in 
Fig. 2 where we plot the piezoresistive anisotropy fac- 
tor (r|| — r_i_)/r|| obtained from the data of Fig.l as a 
function of p — p c . For both bond and site percolation 
models with or without distribution of 2d/t; values, the 
anisotropy factor monotonically goes to zero as p — > p c 
and in the critical region p — p c <C 1 it follows a power 
law behavior: 
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to see from Eqs.(||J|) that in our piezoresistive model 
Py/cx — 1 — (3T 1 1| — r_i_)e, so that the two exponents A and 
A' are equal. 

In summary, we have shown numerically that in 
the critical regime of percolative resistor networks the 
piezoresistive anisotropy has a power law behavior. This 
result suggests a possible experimental route to char- 
acterize transport properties in TFRs or other piezore- 
sistive granular compounds. In fact, as pointed out in 
the introduction, typical TFRs show a power law be- 
havior of the [Conductance G as a function of metallic 
concentrations Hence, we predict that if this behavior is 
due to the closeness to a percolative critical point, then 
the piezoresistive anisotropy factor should follow Eq.(0) 
or equivalently: 
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where t is the critical exponent of the unstrained conduc- 
tance G ~ (p — pc) 1 . We are not aware of published data 
from which it is possible to infer G and (F|| — r_i_)/r|| as 
a function of the metallic concentration in a controlled 
way, so Eq.(||) cannot yet be tested. However, there is 
evidence that the piezoresistive gauge factors anisotropy 
(.Kl — K^/Kl in soma-, commercial TFRs is lower for 
higher resistive samples,u in qualitative accord therefore 
with Eq.(||). We encourage acquisition of additional ex- 
perimental data in TFRs to verify our predictions. 
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This specific power law dependence is made more clear 
in the inset of Fig. 2 where the data are plotted in a log- 
log scale. A fit of the numerical data to Eq.(|7|) leads to 
A = 0.7 ± 0.1 for bond percolation and A = 0.4 ± 0.2 for 
site percolation independently of 2c£/£ within the error 
bars. Although we obtain a critical exponent for bond 
correlation somewhat larger than that for site percola- 
tion, we do not take this difference too seriously since a 
precise determination of A requires much more extensive 
numerical calculations than those presented here. This is 
outside the scope of this communication, which instead 
hints to point out the existence of a critical exponent A 
for the piezoresistive anisotropy. 

The power law behavior of Eq.(Q) can also be inferred 
from the transport properties of anisotrcqaic-bond per- 
colation models studied some time agoOllj In those 
works, percolating networks with random bond conduc- 
tances were defined in such a way that the conducting 
bonds were equal along two directions, for example y 
and z, but different from those along the .-third direc- 
tion, that is x. Topological considerations, til renormaL 
ization group analysis J13 and numerical calculations!!^! 
showed that the quantity o~ y /a x — 1 = a(p) has a power 
law behavior a{p) ~ (p — p c ) x close to the percolation 
thresold. By neglecting the geometrical factors, it is easy 
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